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Abstract 

We investigate the effects of inverse Compton scattering by electrons and positrons 
in the unshocked winds of rotationally-powered binary pulsars. This process can 
scatter low energy target photons to produce gamma rays with energies from MeV 
to TeV. The binary radio pulsars PSR B1259-63 and PSR J0045-73 are both in 
close eccentric orbits around bright main sequence stars which provide a huge den- 
sity of low energy target photons. The inverse Compton scattering process transfers 
momentum from the pulsar wind to the scattered photons, and therefore provides 
a drag which tends to decelerate the pulsar wind. We present detailed calculations 
of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, 
showing that the deceleration of the wind of PSR B1259— 63 due to 'inverse Comp- 
ton drag' is small, but that this process may confine the wind of PSR J0045— 73 
before it attains pressure balance with the outflow of its companion star. We cal- 
culate the spectra and light curves of the resulting inverse Compton emission from 
PSR B1259— 63 and show that if the size of the pulsar wind nebula is compara- 
ble to the binary separation, then the 7-ray emission from the unshocked wind 
may be detectable by atmospheric Cerenkov detectors or by the new generation of 
satellite-borne 7-ray detectors such as INTEGRAL and GLAST. This mechanism 
may therefore provide a direct probe of the freely-expanding regions of pulsar winds, 
previously thought to be invisible. 
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Rotationally-powered pulsars are rapidly rotating neutron stars whose rota- 
tion is slowing at a readily measurable rate. It is generally accepted that the 
rotational energy is carried away by a wind of electrons, positrons, and possi- 
bly ions. Models for the best studied rotationally-powered pulsar - the Crab 
- suggest that the wind expansion is highly relativistic with a bulk Lorentz 
factor of 7w ~ 10^ [Rees & Gunn 1974]. A number of authors have investi- 
gated the inverse Compton scattering of emission from the neutron star itself 
by the electrons and positrons of the pulsar wind both close to the surface 
[Zhang, Qiao & Han 1997; Qiao & Lin 1998; Luo & Protheroe 1998] and at 
about the distance of the hght-cylinder [Bogovalov & Aharonian 1999; Aharo- 
nian & Bogovalov 1999]. At much larger radii, pulsar winds may be confined 
by pressure balance with the surrounding medium, as in the Crab nebula, in 
which case the electrons and positrons in the wind will be accelerated and 
isotropised at the shock which terminates the relativistic wind. The inverse 
Compton emission by electrons and positrons in the shocked pulsar wind has 
been considered, addressing both the scattering of self-produced synchrotron 
photons [De Jager & Harding 1992], of diffuse background target photons [for 
a review see Harding 1996], and of photons from a binary companion [Kirk, 
Ball & Skjaeraasen 1999]. However, until recently, it has generally been thought 
that the intermediate region between the light cylinder and the termination 
shock could not be observed directly. 

In this paper we investigate inverse Compton scattering by the electrons and 
positrons in this unshocked region of pulsar winds. Chernyakova & lUarionov 
[1999] have also considered aspects of this problem, such as the deceleration 
of the pulsar wind when the scatterings can be treated in the Thompson 
limit, and the spectrum of scattered photons when the wind deceleration and 
absorption by photon-photon pair production could be neglected. Here we 
present a more detailed treatment, including deceleration in the Thomson and 
Klein-Nishina regimes and the effects of absorption. We apply our formahsm to 
two binary pulsar systems in which there is a hiigc density of target photons 
from a luminous companion star and show that the pulsar PSR B1259— 63 
should produce detectable gamma-ray emission from the scattering of the 
companion star photons by the unshocked pulsar wind, provided it is not 
confined to a small region around the pulsar. The strong dependence of the 
target photon density on the position with respect to the star, coupled with 
the fact that the target photons may be assumed to all be propagating radially 
from the star, and so are unidirectional at each point where a scattering occurs, 
implies that the scattered emission will have a strong dependence on orbital 
phase. In close binary systems the scattering process can have a significant 
influence on the dynamics of the unshocked wind, and may provide the flrst 
direct probe of the freely-expanding region of a pulsar wind. 

PSR B1259— 63 is the only known galactic radio pulsar which is orbiting a 
main-sequence companion. It is in a close, highly eccentric orbit (e ~ 0.87) 
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around SS2883. a B2c star of radius ~ QRq and luminosity ~ 8.8 x 
10^ L0 [Johnston ct al. 1992, 1994, 1996]. The system is thought to be 1.5 kpc 
from the Earth and at periastron the pulsar is only 23/2* k, 10^^ m from 
its companion. The only other similar binary pulsar system, PSR J0045— 73 
[Kaspi et al. 1994a,b; Bell et al. 1995], is in the Small Magellanic Cloud, and is 
the most distant pulsar known. It is in an orbit of eccentricity 0.81 around an 
even more luminous companion, a Bl V star of radius 6.4i?Q and luminosity 
1.2 X 10"^ L0, and at periastron the pulsar is just 2QRq pa 2 x 10^° m from its 
companion. 

The spin-down luminosity of PSR B1259— 63 is significant, Lp = 8.3 x 10^^ W, 
and it is likely that this drives a relatively strong wind. The wind is most 
probably confined by pressure balance with the strong outfiow from the pul- 
sar's Be-star companion [Melatos, Johnston & Melrose 1995], in which case 
the pulsar wind slows to a subsonic outflow at a termination shock. The al- 
ready relativistic electrons and positrons from the pulsar wind are isotropised 
and accelerated to even higher energies at this shock. The unpulsed X-rays 
observed throughout the binary orbit of this system [Hirayama at al. 1996] 
are thought to result from synchrotron emission from electrons and positrons 
in the shocked region of the pulsar wind [Tavani & Arons 1997]. Kirk et al. 
[1999] have shown that inverse Compton scattering of photons from the Be 
star, by the electrons and positrons of the shocked pulsar wind, is likely to pro- 
duce a flux of hard 7-rays that would be detectable by atmospheric Cerenkov 
telescopes. 

If just 10^'^ of the spin down luminosity of PSR B1259— 63 is converted into 
7-rays by inverse Compton scattering in the unshocked wind, the resulting flux 
would be detectable at energies somewhere in the MeV-TeV range depending 
on the Lorentz factor of the pulsar wind. The lower spin down luminosity and 
greater distance to PSR J0045— 73 means that it should be unobservable in 
hard 7-rays. Nevertheless the binary separation of this system is even smaller 
than that of PSR B1259— 63 and the companion star is also more luminous, so 
inverse Compton scattering by the unshocked wind may affect the dynamics 
of this system. 

In the following sections we present equations describing the deceleration of 
the pulsar wind by 'inverse Compton drag' - the energy losses of the unshocked 
electrons and positrons due to inverse Compton scattering of the companion- 
star photons, and the intensity and spectrum of the energetic scattered pho- 
tons. We show that the inverse Compton losses cannot contain the wind of 
PSR B1259— 63 before it attains the radius at which pressure balance with 
the Be-star outflow occurs. These losses may signiflcantly decelerate the wind 
of PSR J0045— 73 at times near periastron within a range of angles about the 
line joining the pulsar and its companion. We present examples of the emit- 
ted inverse-Compton spectrum and of the light curves from PSR B1259— 63, 
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showing that the emission from the unshockcd wind should be detectable as 
gamma-rays. The signal from the unshocked wind is quite different from that 
from the shocked wind, so that if a variable gamma-ray signal is detected from 
this system, these two possible source mechanisms should be readily distin- 
guishable. 



1 Dynamics of the pulsar wind 

We assume that a relativistic, radially directed wind containing electrons, 
positrons and an admixture of ions emerges from the pulsar. This wind prop- 
agates in an ambient radiation field and the particles in the wind interact 
incoherently with the ambient photons. The most important interaction is 
Compton scattering by the electrons and positrons. In this section we outline 
the essential steps in the derivation of the equations for the dynamics of a pul- 
sar wind subject to inverse Compton scattering. Further details are presented 
in Appendix A. 

Observations and models of the wind of the Crab pulsar suggest that it emerges 
nearly tangcntially at around the light cylinder radius and that the Poynting 
fiux dominates the outfiow at such radii [Rees & Gunn 1974]. However, by 
the time the wind interacts with the surrounding nebula the importance of 
the magnetic field is greatly diminished and the outgoing fiux is dominated 
by the particle momentum. Models suggest that the wind ultimately attains 
a bulk Lorentz factor of order 7^ = 10^ [Kennel & Coroniti 1984; Emmering 
& Chevalier 1987]. While this picture is well accepted, the physics responsible 
for the transformation of the wind is poorly understood. Nevertheless, it is 
likely that the wind evolves relatively rapidly into a cold, relativistic, particle 
dominated, radial outfiow. We will assume that at the radii of interest here, 
which are well beyond the light cylinder, the electrons and positrons of the 
pulsar wind may be taken to be monoenergetic, with Lorentz factor equal to 
that of the bulk fiow, 7w 

— * — * 

We denote by n^{e, R, Q) dQde the differential number density of target pho- 
tons moving within the solid angle dQ of the unit vector Q at position R with 
respect to the star, with energy between E = emeC^ and (e + de)mcC^. We 
investigate the deceleration of the wind in the highly anisotropic radiation 
field comprising the light from a companion star which is treated as a point 
source. The target photons are treated as monochromatic. 
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An individual electron of Lorentz factor 7 scatters target photons at an average 
rate 



Nic = c dQ dcin — n^(ein,R,fi)aKN(e( 
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[Jones 1965], where crKN(a;) is the Klein-Nishina cross section. (Note that when 
used alone the term 'electron' may be taken to refer to both electrons and 
positrons in the pulsar wind.) Unprimed quantities are measured in the ob- 
server's frame and primed quantities are measured in the rest frame of the in- 
coming electron. Energies are expressed in units of rrieC^ so that i?in = einW^eC^ 
is the energy of the target photon in the observer's frame, = e[^me(? the 
energy of the target photon seen in the electron rest frame. The process of 
triplet pair production, e7 3e, can be ignored since the maximum value of 
7ein which we consider is of the order 100 [Mastichiadis 1991]. 

The average rate of change of energy of the photon field is given by an integral 
similar to that in Eq. (0), but weighted by the factor eout — Cin, where i^out = 
Cout^^eC^ is the energy of the scattered photon given by 



47(1-/3 cos Ct) 
1 -u ^' n _ nr^^ni \ 



1 + 4(1- COS 

with 6'^^^^ the angle between the outgoing photon and the velocity vector of the 
observer's frame and 9'^^ the scattering angle, both measured in the electron 
rest frame. The rate of change of radial momentum is calculated in a similar 
fashion, using —eout cos 6'out + Cin cos 6^^ as the weighting function, with 6*10 and 
^out the angles between the incoming and outgoing photons and the vector 
opposite the incoming electron velocity in the observer's frame. 

As a result of Compton scattering the electrons and positrons of the pulsar 
wind lose both energy and radial momentum, but gain a small amount of 
transverse momentum, so that the gas is both decelerated and heated. In the 
case of an ultra-relativistic wind the heating is small and may be neglected 
(see Appendix A). The deceleration of the wind due to inverse Compton drag 
is then given by the equation 

^ = / d^de;,^ (l-l^-'J^\ n^(e,^, R, moMn) (3) 

dr J Cin V 7ein 7 / 

in which = (1 — Iw'^Y^'^j '^t is the Thompson cross section, and Fioss is a 
function defined by Jones [1965] 

-2x(10x^ - 51x=^ - 93x2 _ _ 



3x^(1 + 2x)3 
(x2 - 2x - 3) ln(2x + 1) 



(4) 
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[see also Kirk et al. 1999]. The parameter t] is the fraction of the wind lumi- 
nosity carried by electrons and positrons: 



^ 1 + Mn-J{m^ne) ^^'^ 

with M the ion mass and me the electron mass, and n\ and rie the ion and 
electron number densities at any given radius. 



pulsar 




Earth 



Companion star 

Fig. 1. Sketch of the binary system defining angles and distances used in the text. 

When the target photons originate from a companion star, the target field can 
be written in the form 

n^(ei„ R, VI) = N{R)5{t,^ - t,)5{Vl - R) (6) 

where N{R) = L^/ {AT^B^chvo) is the photon density, eo = huQ/nieC^ is the 
dimensionless energy of the companion star photons, R = R/i? is a unit 
vector, and is the stellar luminosity. We denote the binary separation of the 
pulsar and its companion star by D, and introduce the dimensionless distance 
f = r/D. Inserting Eq. @ into @, making the approximation /^^ ~ 1, we 
find 

d7w _ SrjNoaTD e'J sin^ ^ / eg \ , 

dr Ssm^^p eo V 7w4 7w/ 

where Nq = N{D), = 7weo(l — /^wCos-?/') and ip is the angle between the 
vectors R and r. The geometry is shown in Fig. |l]. Note that the angle ip is a. 
function of position f: 

f — cos 9p 

cos ip ^ 



1 + — 2r cos 6p 



Equation (0) implies that 7„ > when e[„ < eo/(7w — eo) meaning that the 
target photons actually accelerate electrons traveling almost directly away 
from the companion star for which cosip ~ 1. The condition for this to occur 
is 

cosV^>^^^= (l -) ; (9) 

^1 - 1/7^ V 7w(7w-eo)y 



6 



for 7„ = 10^ and eo = 10~^ this is equivalent to |'?/^|<10~^, which is satis- 
fied only when the electrons are moving essentially radially away from the 
companion star. 



2 An approximate treatment 



Before applying the results of §1 and deriving the spectrum of the scattered ra- 
diation, we consider a greatly simplified treatment obtained by approximating 
the scattering cross-section by the Thompson cross-section and assuming that 

the energy distribution of the scattered photons is a delta function. In this 
approximation an electron of Lorcntz factor 7 emits monochromatic photons 
of energy eout = 47^ein/3 and the reaction rate is simply Na^c. The equation 
governing the deceleration of the wind is then 

^^—Vlle^nNaT. (10) 

If the density of target photons remains roughly constant on the length scale 
upon which the wind is decelerated, it follows that 

rdecel ~ [?77w(0)einA^(7T]"^ (11) 

= 7.8xl0M meters, 12 

where the scahng of the energy density of the target photons has been chosen 
to be close to that of PSR B1259— 63 at periastron. 

We next define a dimensionless spectral energy distribution J-'{r), which de- 
scribes how the radiation builds up along the line of sight within the pulsar 
wind. At the outer boundary r = rmax of the source, the final, observable vahic 
is attained: JF(rmax) = EFE'^i^d'^ / where Fe is the energy flux density at 
the Earth of photons of energy and d is the distance to the pulsar. The 
equation for the generation of the spectral energy density is then 

^ = Na^el^,5 {e,^, - ^ile-ji) . (13) 

Under the same assumption that the density of target photons is roughly 
constant over the distance raeceb the emitted spectral energy distribution is 

(14) 

extending from a maximum at eout = e^Jut = 47^(0)ein/3 down to a minimum 
which is a function of the physical extent of the wind. 
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The approximations given by Eqs. (JT^) and ([T^) can be applied to PSR B1259— 63 
and PSR J0045— 73 at periastron using the appropriate physical parameters 
which are summarized in Table 1. The energy density in the target radia- 
tion field near the surface of PSR B1259— 63 at the time of periastron is 
roughly 6 x 10^^ eVcm~^, which implies rdecei ~ 1-3 x 10^*^ m ~ 0.1 Dr-, where 
Dr = 9.6 X 10^° m is the separation of the stars at periastron. Taken at face 
value this estimate suggests that the inverse Compton drag may be so effi- 
cient as to decelerate the pulsar wind completely before it attains the radius 
at which it would otherwise reach pressure balance with the Be star outfiow. 
Furthermore this estimate implies that there should be a substantial energy 
fiux at the Earth (~ lO^eV cm s ) concentrated in TeV photons. This cor- 
responds to a photon fiux of ~ 10"^ cm^^ s~^, roughly 10^ times higher than 
the sensitivity threshold of the imaging Cerenkov detector operated by the 
CANGAROO collaboration, and clearly confiicts with the marginal detection 
found by that experiment [Sako et al. 1997]. 

For PSR J0045— 73 the energy density in the target radiation field at periastron 
is even higher, [/rad ~ 2 x 10^^ eV cm~^, whence raecei ~ 1-3 x 10^° m ~ 0.02 Dr 
with D-r = 1.8 X 10^° m for this system. However, the much greater distance 
of this source makes it uninteresting for gamma-ray observations. 

These simple estimates indicate that inverse Compton drag on the unshocked 
winds of these binary pulsars may be important. However the Klein-Nishina 
corrections to the scattering cross section are significant in these systems, as 
are the effects of the changing angle ip and the changing density of target 
photons encountered by the scattering electrons as they propagate radially 
from the pulsar. We include these effects in the next sections. 
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Tabic 1 

Parameters for the PSR B1259— 63 and PSR J0045— 73 systems. References are: 
1-Johnston et al. [1996]; 2- Johnston et al. [1994]; 3-Underhill & Doazan [1982]; 
4-Kaspi et al. [1994a]; 5-Kaspi et al. [1994b]; 6-Bell et al. [1995]. 



Parameter 
Pulsar 

period P 
period derivative P 
surface magnetic field 
spin down luminosity Lp 
Companion star 
spectral type 

effective temperature T^s 
radius i?* 

luminosity 



PSR B1259-63 



Value 

47.762 ms 
2.279 X 10"^^ 

3.3 X lO'^ T 
8.3 X 10^8 W 

B2 

2.28 X 10^ K 

6.Oi?0 
(4.2 X 10^ m) 
8.8 X lO^L© 



(3.3 X 10^0 W) 

effective photon energy eo = 2.7kB,TeS / {mc?) 

10-^ 
lOM© 
(2 X 10^1 kg) 

e 0.87 
Dt 23i?* 

(9.6 X 10^° m) 
Da 2>ZIR* 
(1.4 X 10^2 

i 35° 
d 1.5 kpc 



mass 

System 

eccentricity 
periastron separation 

apastron separation 

orbital inclination 
distance 



Ref. 

1 
1 
1 
1 

2 
3 
3 



PSR J0045-73 
Value Ref. 

0.926 s 4 

4.486 X 10^^^ 5 

2.1 X 10^ T 4 
2.2 X 10^5 W 

Bl V 6 

2.4 X 10^ K 6 

6.4ii:0 6 

(4.5 X 10^ m) 

1.2 X IQ'^Lq 6 
(4.6 X 10^0 W) 

10-5 

8.8M0 6 
(1.8 X 10^1 kg) 

0.81 4 

(1.8 X 10^° m) 

38i?* 
(1.7 X 10^1 m) 

44° 6 
57.5 kpc 4 
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3 Inverse Compton emission from the unshocked wind 



The basic approximation made in computing the inverse Compton emission of 
a relativistic particle is that the scattered photons are emitted in the direction 
of motion of the particle. Thus in the case of scattering by a radial pulsar wind, 
an observer sees photons produced at all points along the radius vector from 
the pulsar to the point at which the wind is either decelerated by radiation 
drag, or shocked, or simply runs out of target photons. In addition, scattered 
photons of sufficiently high energy may be absorbed by pair production on the 
same population of target photons via the process 77 e~^e~. The probability 
for pair production, like the scattering rate, is a function of position along the 
radius vector. In this section we derive expressions for the inverse Compton 
emission from the unshocked wind, including the effects of pair-production 
absorption. Further details are presented in Appendix B. 

Observations of the Crab nebula suggest that in that case at least, the pulsar 
wind has an energy-dependent anisotropy [Aschenbach & Brinkmann 1975; 
Hester et al. 1995]. However, for simplicity we assume that the pulsar wind is 
independent of latitude and longitude. The distribution function of monoen- 
ergetic, radially directed electrons in the wind may then be written as 



where iVe(r) = rjL^/ [47rr^c/3„7w(0)meC^], and L„ is the luminosity of the 
pulsar wind, which is taken to be equal to the spin-down luminosity. It is 
relatively straightforward to incorporate a more complicated model, provided 
the flow remains radial. 

The specific intensity of radiation /(eout, 1^) is governed by the radiative trans- 
port equation 



dr \AtcJ 7w(0)mcC^ deoutdt dr' 

where diV^/deoutdt is the differential rate of emission of inverse Compton scat- 
tered photons by a single electron, and r is the optical depth between the 
observer and the point r due to pair production. Since r is measured from 
the pulsar outwards, it follows that dr/dr < 0. More details are presented in 
Appendix B. 

When the target radiation originates from a companion star, as specified in 
Eq. (^, the photon production rate can be written as 



ne(7, r, n) = N,ir)5 [7 - 7w(r)] 6(^1 - r) 



(15) 




dK 
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N{R)aTC 7Y(eout, cosip), 



(17) 



dtde, 



out 
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where 7i(eout, cosip) is related to the quantity dA^-^/dtdeout[deindx'] derived by 
Ho & Epstein [1989] and is defined in Appendix B. 

In terms of the dimensionless radius f = r/D, the dimensionless spectral en- 
ergy distribution J^(r) introduced in the previous section (jF(oo) = And'^EFE/L^ 
where Fe is the flux and d the distance to the object) builds up in the pulsar 
wind according to: 



where Bp is the angle between the line joining the stars and the line of sight. 
The absorption optical depth obeys the equation [Kirk et al. 1999] 



where cr^'y{z) is the cross section for pair production in units of ctt (see Ap- 



pendix B), which is to be evaluated with argument e = A/eouteo(l — cosip) /2. 



The computation of the inverse Compton luminosity of the pulsar wind in the 
radiation field of its companion star therefore reduces to the integration of 
equations (^ and (0), with the optical depth gradient given by ([l^). 



4 Inverse Compton drag on the winds of 
PSR B1259-63 & PSR J0045-73 

Equation (|^) can be numerically integrated to obtain the Lorentz factor of 
the radial wind as a function of r. Figure |^ shows the results for four different 
values of 9p, for parameters appropriate to PSR B1259— 63 and PSR J0045— 73 
as given in Table [l|, assuming that the wind momentum is carried entirely by 
electrons and positrons, i.e. = 1. When = the electrons are traveling 
straight towards the companion star and inverse Compton drag is at its most 
effective because the collisions are head on and the density of target photons 
increases as (1 — f)~^. For large values of 6p the drag is relatively ineffective 
because the electron-photon collisions are far from head on, and furthermore, 
the target photon density drops as the electrons propagate away from the 
companion star. 

The curves for = shown in Fig. ^ can be compared to the simple estimates 
presented in §0. Figure ^ shows that in the head-on case, inverse Compton 
drag stops electrons in the wind of PSR B1259— 63 at f 0.9 compare to the 
simple estimate rdecei ~ O.lDr. On the other hand, Fig. 0b shows that head-on 




(18) 




(19) 
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Fig. 2. Lorentz factor 7w(^) of a radial wind with 7w(0) = 10^ and rj = 1, subject 
to inverse Compton drag from the companion star photons, for PSR B1259— 63 (a) 
and PSR J0045— 73 (b). The pulsar-companion star separation is that at periastron. 
The curves correspond to different values of 9p (see Fig. as labelled. 

scatterings stop the electrons in the wind of PSR J0045— 73 at f ^ 0.7, com- 
pared to the simple estimate rdecei ~ 0.02/^,-. These differences arise primarily 
because the simple estimates ignored the Klein-Nishina effects — which re- 
duce the scattering cross section from the Thompson value. The increasing 
density of target photons encountered by the wind electrons as they approach 
the companion star offsets the Klein-Nishina correction to some degree, but 
the net effect is that the deceleration length is much larger than suggested by 
the simple estimate. 

Equation (J^) can also be used to determine the position of the surface at which 
the radial wind loses a given fraction of its energy as a result of inverse Comp- 
ton scattering. Such a surface is rotationally symmetric about the line joining 
the two stars. Cross sections of the surfaces where 7w('^) has fallen to 90%, 
75% and 50% of its initial value are shown in Fig. ^ for PSR B1259— 63 and 
PSR J0045— 73 at periastron, for three different values of the initial Lorentz 
factor 7„(0). Unless the interaction between the stellar wind and the pul- 
sar wind intervenes, we expect the pulsar wind to terminate roughly at the 
position where the wind electrons have lost half of their initial energy. 

The shapes of the cross-sections in all the panels of Fig. ^ clearly show the 
strong dependence of the wind deceleration on the angle 9p. The surfaces of 
equal loss are all closest to the pulsar along the line joining the pulsar to the 
Be star, and rapidly become more distant from the pulsar as 6p increases. 
This occurs because while head-on collisions are effective at decelerating the 
electrons, the deceleration drops off very rapidly as 9^ increases, and the pho- 
ton density also drops off very rapidly with distance from the pulsar at even 
moderate values of 9p. 

The dependence on 7w(0) as depicted in Fig. ^ is particularly striking. The 
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Fig. 3. Cross sections through the surfaces where 7(r) has dropped to 90%, 75% and 
50% of the initial value 7w(0) - as labelled, if inverse Compton drag is the dominant 
wind deceleration mechanism. The wind, orbit and target photon parameters are as 
for PSR B1259— 63 (left) and for PSR J0045— 73 (right) at periastron and are given 
in Table The pulsar is at the origin and is marked by a plus and the companion 
star is marked by an asterisk. The two top panels (a) & (b) show the results when 
7w(0) = 10^, the middle two (c) & (d) are for 7w(0) = 10^, and the lower two panels 
(e) & (f) are calculated for 7„(0) = 10*^. 
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surfaces of equal loss are much closer to the pulsar, for both PSR B1259— 63 
and PSR J0045-73, when 7^(0) = 10^ than for either 7^(0) = 10^ or 7„(0) = 
10^. This in an indication that the deceleration is most effective when the 
system is 'tuned' such that 760 ~ 1. This occurs because the energy loss rate 
is made up of two parts, the scattering rate and the average energy loss per 
collision. The former is constant in the Thomson regime, but falls off with 
increasing energy in the Klein-Nishina regime. The average energy loss per 
collision divided by the energy of the scattering particle increases proportional 
to 7 in the Thomson regime, and tends to a constant in the Klein-Nishina 
regime. The combined result of these effects is that there is a peak in the 
energy loss rate near the transition between the Thomson and Klein-Nishina 
regimes. 

Finally, Fig. ^ shows the difference between the inverse Compton deceleration 
the winds of PSR B1259-63 and of PSR J0045-73. In PSR B1259-63 the 
density of photons from the companion star is not large enough to substantially 
decelerate the pulsar wind before it has travelled a substantial fraction of D^, 
even along the line joining the stars where the collisions are precisely head 
on. The surfaces depicted, which all have their apexes at f>0.25, therefore 
wrap around the Be star rather than the pulsar. These results suggest that in 
PSR B1259— 63 the inverse Compton drag is unlikely to terminate the pulsar 
wind before it attains pressure balance with the Be star outflow, given that 
this balance is expected to occur at f <0.5 [Melatos et al. 1995]. It is therefore 
reasonable to assume, with Kirk et al. [1999], that a termination shock forms 
in the pulsar wind as a result of pressure balance with the Be-star outflow. 
Nevertheless, in the case of PSR B1259— 63 the conversion of even a small 
fraction of the wind energy could result in an observable flux. From Fig. |^ it is 
clear that if the wind terminates at a radius tt which is a significant fraction 
of the binary separation, it may become observable. For 7w(0) ~ 10^ and 
fx ~ 0.5 the figure indicates that at periastron the unshocked wind will lose 
approximately 25% of its momentum due to inverse Compton scattering along 
the line joining the stars, resulting in a substantial gamma-ray luminosity. 

For PSR J0045— 73 the situation is even more extreme. The companion star 
photon density at the pulsar is so large at periastron - because the orbital 
separation is so small - that inverse Compton drag is much more effective at 
decelerating the pulsar wind. In this system many of the surfaces of equal loss 
depicted have their apex very close to the pulsar and therefore wrap around 
the pulsar rather than the companion star. It is possible that inverse Compton 
drag may stop the wind of PSR J0045— 73 before it attains pressure balance 
with the companion star outflow, at least at angles close to the line connecting 
the stars, and at times close to periastron. 
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5 Emitted gamma rays 



5. 1 Spectra 

This section details the spectrum of inverse Compton emission produced via 
scattering from the unshocked wind. We assume that the wind momentum is 
carried entirely by electrons and positrons [r] = 1), that the only deceleration 
is that from the scattering, and that the wind expands freely without being 
terminated via pressure balance with the companion star outflow {jt ^ !)• 
The spectra obtained by integrating equations (|^ and (|T8|) with (|19|), for a 
range of wind Lorentz factors 7w(0) and at three different angles 6'p with other 
parameters relevant to PSR B1259— 63 at periastron, are presented in Fig. ^. 
Note that the spectra show the energy flux versus the photon energy. The 
emitted spectra depend dramatically on both the angle 6^ and on 7w(0). 

For 6'p = 5° (a value not, however, observable from the Earth) the wind elec- 
trons are initially directed almost along the line connecting the pulsar and 
its companion. Inverse Compton scattering is thus intrinsically efficient as the 
electron-photon collisions are close to head on, and the electrons in the wind 
encounter a very large target photon density as they pass close to the com- 
panion star. Integration over the scatterings occurring along such a radius 
therefore produces a relatively broad spectrum of emitted photons because 
of the wide range of Lorentz factors of the scattering electrons. Furthermore, 
the scattered emission is relatively bright because the wind loses a substantial 
fraction of its momentum to the scattered photons. 

For 9p = 125° the wind electrons are initially directed more or less away from 
the companion star, so the scatterings are very far from head on and the target 
photon density decreases very rapidly as the electrons propagate away from 
the pulsar (and hence from the companion star). The emitted spectrum in this 
case is close to monochromatic because the electron deceleration is very small 
so the scattering particles are essentially monoenergetic. Furthermore, since 
only a small fraction of the wind momentum is transferred to the scattered 
photons, the flux of the scattered photons is commensurately small. 
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Fig. 4. The SED (spectral energy distribution) EFe with Fe the energy flux at the 
distance of the Earth of inverse Compton scattered photons as a function of photon 
energy - at angles 9p = 5° (left), 65° (middle), and 125° (right). The top row of 
spectra are calculated for a pulsar wind with an initial Lorentz factor 7w(0) = 10^. 
Each subsequent set of three spectra down the page has 7w(0) increased by a factor 
of 10, so that the bottom row represents the emission for 7„(0) = lO'^. The other 
parameters are the same for all panels, and are appropriate for PSR B1259— 63 at 
periastron. The dashed curves show the spectra calculated by ignoring the absorp- 
tion due to pair production, i.e. by artificially setting dr/dr = in Eq. ([l8|). The 
solid curves include the effects of the pair production optical depth. Note that the 
energy scales are the same for each set of three spectra for a given wind Lorentz 
factor, but the scale shifts up a decade with each factor of ten increase in 7w(0). 
The energy flux scales are the same for all panels. 

The dependence of the emitted spectra on 7w(0) is also striking. The energy of 
the scattered photons increases strongly as the wind Lorentz factor increases; 
e.g. for 7w(0) = 10'^ the peak of the emitted spectrum occurs just above 10 MeV 
while for 7w(0) = 10^ the peak is at ~ 1000 MeV. The energy of the scattered 
photons therefore scales as ~ 7^ in this regime, as expected since for these 
Lorentz factors 769 -C 1 and the scattering is in the Thompson regime. As the 
wind Lorentz factor is increased further the peak of the spectrum of scattered 
photons continues to move up in energy, but the increase is slower than 7^ due 
to the increasing importance of the Klein-Nishina effects. The broadening of 
the spectrum as 7w(0) increases to around 10^ and the subsequent narrowing as 
7 increases further is dramatically illustrated by the left-hand panels for 6p = 
5°. This is another example of the effects of tuning the system to maximize the 
scattering efficiency, resulting in the greatest deceleration and correspondingly 
broad energy distribution of the scattering particles, at 769 ~ 1. 

The shape of each individual spectrum is also indicative of the degree of de- 
celeration effected by the scattering. All the spectra cut off abruptly at high 
energies reflecting the maximum energy e^^t of the scattered photons. The 
spectrum calculated for 7^ = 10^ and 9p = 5° rises relatively steeply at the 
lowest energies shown, then turns over and rises more slowly for about 1.5 
decades in energy (above 10'^ MeV) before dropping steeply towards the upper 
cutoff. The steeply-rising portion of the spectrum at low energies is increasing 
as e^ut corresponding to a flux density which decreases as 6^,^/^ - the usual 
spectrum for inverse Compton emission from a single electron in the classical 
Thompson regime. This portion of the emitted spectrum is dominated by the 
electrons at the minimum Lorentz factor attained after deceleration of the 
wind. The flatter portion of the spectrum corresponds to the emission from 
the decelerating portion of the wind — integrated over a relatively broad range 
of 7^; it is only evident when both 7w(0) and 9p have values which lead to a 
substantial deceleration. When this is not the case the spectra take one of two 

1 /2 4- 

forms: if 760 -C 1 the SED increases as e^ut right up to the cutoff at e'^^^; if 
760 ^ 1 the Klein-Nishina effects are important and the SED rises much more 
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steeply to a peak near e^nt and the emission is essentially monochromatic. 

Finally, Fig. ^ shows the effects of absorption due to pair production. For 
7w(0) <10^ the scattered photons are below the threshold energy required for 
pair production on the target photons even when 6p is small, and so there is no 
absorption. The absorption due to pair production is most significant in the 
spectrum plotted for 7w(0) = 10^ and 9p = 5° where it removes a substantial 
part of the emission at energies up to a decade below the peak. For smaller 
7w the scattered photons are below the threshold for absorption, whereas at 
higher values the resonant character of the cross section becomes apparent. 



5.2 Light curves 

The inclination angle of the orbit of PSR B1259— 63 is i = 35° so a telescope 
on the Earth will sample angles 6p between 90 — i = 55° and 90 + i = 125° 
over the binary period. These extremes are the two larger values of 6p for 
which the spectra are plotted in Fig. ^. The increase from the minimum value 
of 9p to the maximum occurs over just 65 days beginning on day —6, just 
before periastron. The subsequent decrease occurs far more slowly over the 
remaining 1170 days of the orbital period. This variation, together with the 
variation in the target photon density due to the changing binary separation 
during the orbit, results in a profound variation in the emission of inverse 
Compton scattered radiation from the binary system. 

Fig. |] shows the orbital variation of the integrated flux density, / FsdE, 
expected from PSR B1259-63 for 7^(0) = 10^ 10^ and 10^. The light curve 
for 7w(0) = 10^ is very similar to that for 7w(0) = 10^ so it has been omitted 
for clarity. The total variation in the integrated flux is very large; for 7„ = 10^ 
the peak is two orders of magnitude above the minimum. The light curves are 
also very strongly peaked and clearly asymmetric, with the increase towards 
the maximum integrated flux occurring somewhat more gradually than the 
decline to the minimum. The asymmetry, peak to minimum ratio, and peaked 
nature of the light curves are minimized when 7„(0)eo ~ 1, but even then the 
peak to minimum ratio is more than a factor of ten. These features are all a 
direct reflection of the orbital variation of the binary separation and the angle 
Op] the asymmetry is a straightforward manifestation of the asymmetry in the 
variation of 6p about periastron. 

The peak of the spectrum of the scattered photons also varies with orbital 
phase. If this variation were to cross the threshold of an observing instrument 
the detected light curve could be somewhat different from those shown in 
Fig. ^ The energy spectra at the extreme values of 6p sampled during one 
binary orbit are shown in the second and third columns of Fig. ^. The variation 
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Fig. 5. Light curves showing the integrated energy flux at the Earth from 
PSR B1259— 63. The different curves correspond to different initial wind speeds: 
the solid curve is for 7w(0) = lO'*, the dashed curve is for 7w(0) = 10^, and the 
dot-dashed curve is for 7w(0) = 10''. The left panel shows the variation over the 
entire orbit and the right panel shows just 200 days centered on periastron. 

in the peak energy is largest when 7^(0)60 ^1. For 7w(0) = 10^ the peak is 
at ~ 1000 MeV just prior to periastron (when 6p attains its minimum value), 
and it moves down to ~ 200 MeV by about day +60 (corresponding to the 
maximum of ^p). For 7w(0) = 10^ the peak moves from ~ 30 GeV to ~ 20 GeV, 
while for 7^(0) = 10*^ and 7^(0) = 10^ the peak is at ~ 300 GeV and ~ 5TeV 
respectively, and moves very little. 

The integrated flux shown in Fig. |^ does not include the effects of absorp- 
tion due to pair production. Calculations including the pair production optical 
depth indicate that its effect on the integrated flux density from PSR B1259— 63 
is small - reducing the integrated flux density by ~ 10% or less - and that the 
reduction is relatively insensitive to the orbital phase. It is therefore unimpor- 
tant for observations of the inverse Compton emission. 
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6 Discussion 



The models for the inverse Compton emission from the freely-expanding wind 
of PSR B1259— 63 presented in §5 suggest that the emission is hkely to be 
detectable at energies somewhere in the GeV to TeV range, depending on the 
wind parameters. Given the present lack of operational 7-ray telescopes, the 
best opportunity for detecting the emission is in the TeV energy range. 

In calculating the inverse Compton emission spectra we optimistically assume 
that none of the wind momentum is carried by ions and that the wind is not 
terminated via pressure balance with the companion star outflow. The results 
indicate that if the initial Lorentz factor of the wind is in the range 10^—10'' the 
inverse Compton emission from the unshocked wind at epochs near periastron 
could be as much as a factor of 100 — 1000 above the sensitivity threshold of the 
new CANGAROO II imaging Cerenkov detector [Yoshikoshi et al. 1999] which 
should be operational before the end of 1999. If the pulsar wind is slower, the 
inverse Compton emission should be detectable by the GLAST observatory 
at energies in the range 20MeV-300GeV if 7^(0) ~ 10^ - 10^ or by the 
INTEGRAL project at energies up to lOMeV if 7w(0) ~ 10^. The sensitivity 
of each of these planned satellite-borne observatories is easily sufficient to 
detect the model emission. However, they will not operate until well after the 
next periastron. 

If either of the fundamental assumptions regarding the pulsar wind is not 
applicable to the PSR B1259— 63 system, the gamma-ray emission from the 
unshocked wind will be lower than suggested by the models detailed in §5. 
If the momentum flux of the pulsar wind is sufficiently small that it is dom- 
inated by the Be-star outflow, then the pulsar wind will terminate close to 
the pulsar at a shock which is wrapped around the pulsar. The mismatch be- 
tween the shape of this termination shock and the loss surfaces due to inverse 
Compton scattering could then lead to integrated-flux light curves that are 
quite different to those presented in Fig. |^. Furthermore, the inverse Compton 
deceleration of the pulsar wind decreases its ram pressure, reducing the termi- 
nation shock radius, but this effect will be small at values of Op which can be 
observed. The nature of the gamma-ray light curves will depend on both the 
stand-off distance and the details of the shape and orientation of the pulsar 
wind termination shock. 

Observations of changes in the pulsed emission from PSR B1259— 63 [John- 
ston et al. 1996] and observations and interpretation of the unpulsed radio 
emission seen around periastron [Johnston et al. 1999; Ball et al. 1999] indi- 
cate that the pulsar and its wind interact with a disk surrounding the Be star 
between about day —25 and +25. There are few observations of X-ray and 
7-ray emission from this system at epochs outside this range. Tavani & Arons 
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[1997] have addressed the problem of fitting the X-ray observations close to 
periastron with the inverse Compton emission from the shocked pulsar wind. 
Their analysis indicates that the details of the flow pattern, the competition 
between adiabatic and radiative losses, and strong shielding (or shadowing) 
of the target photons by the pulsar disk, are all important near the epoch of 
periastron. These details are not considered in our work here, and the light 
curves of Fig. |^ can therefore not be expected to be accurate at epochs close 
to periastron. Nevertheless the upper limits from EGRET measurements [Ta- 
vani et al. 1996] just prior to the 1994 periastron place some constraint on the 
pulsar wind parameters. In particular, the inverse Compton spectrum from 
an unbounded lepton wind with an initial Lorentz factor of 7w(0) = 10^ (as 
shown in Fig. ^ would exceed the EGRET upper limits which imply that 
the energy flux of photons at energies of 30 MeV to 10 GeV is less than about 
2 X lO-^MeVcm-^s-i. 

The periastron to apastron ratio of inverse Compton emission from an un- 
bounded lepton wind from PSR B1259— 63 as shown in Fig. ^ is much larger 
than the flux ratio from the shock-accelerated electrons when the wind termi- 
nates close to the pulsar [Kirk et al. 1999]. Observations of the high-energy 
light curve may therefore enable the deconvolution of the two contributions 
to the 7-ray emission. In particular, observations aimed at detecting an asym- 
metric light curve at energies in the MeV-TeV range, or at least placing firm 
upper limits on the high energy emission from this system, should provide 
valuable information on the properties of the unshocked pulsar wind. 
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Appendix A: Energy and momentum balance 



Assuming the magnetic field to be weak, the equations of conservation of 
energy and radial momentum flux in the wind take on the forms: 



1 d 

J.2 ^j. 



r^(£) + P),8.,T^]=C» s (20) 



ii:p[(B + P)/jj7^ + P]} = C' (21) 

[e.g., Landau & Lifshitz 1959] Here, E and P denote the proper gas energy 
density (including rest-mass) and pressure, respectively. The radial velocity of 
the gas is c/3w and the associated Lorentz factor is 7^. The components C° and 
C^' of the four-force density on the right-hand sides of Eqs. ( pOD and (^) are 



proportional to the proper number density Ue of electrons and positrons taken 
together and to the average rates of energy and radial momentum change of a 
single electron due to Compton scattering: (de/dr) and (dpr/dr) where r is 
the proper time, measured in the rest frame of the wind. 

The average rate of change of energy of the photon field is given by an integral 
similar to that in equation ([l|) with the additional weight factor eout — Cin- 
Equating the energy loss rate by the electron to the gain by the photon field 
leads, after a straightforward calculation, to the expression 

o ^in \ 7^in 7 / 

r^^(ein,R,^i)Floss(e:J . (22) 




The rate of change of radial momentum is calculated in a similar fashion, using 
the quantity —eout cos 6'out + ^in cos 9^^ as the weighting function, which gives 




3 /■ ,^ , e' ' 



-ctt^tt-cC J ddde 



111 

in 




8p J tin V lin 

n^(ein,R,ll)Fioss(e;n) • (23) 



The amount of heating is given by the projection of the four-force density 
onto the velocity. Inspection of Eqs. (|22|) and (^) reveals that only the fourth 
term e(n/7^ in parentheses in Eq. (P3|) survives the projection. In the case 



of an ultra-relativistic wind this term is small (~ ^m/l) compared to the 
leading ones in the expression for the work done decelerating the electron, 
and we neglect it in the following. In this approximation, we set P = and 
E/c^ = p„, the proper rest-mass density, and identify the Lorentz factor of the 
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wind with that of each individual electron. As a result, one of the equations 
([20| ) and (|2T|) is superfluous. Keeping the energy equation (|20|) , and combining 
it with the equation of mass conservation, which reads 

4:r(^V/5w^w)=o (24) 

one finds 

d7w _ V /de 
dr rUeC^P^'jy, \dr 

T] I de 

rrieC^f^w \ dt 

Equation follows directly through substitution of (^) into (pSj). 
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Appendix B: Radiation transfer of inverse Compton emission 



The differential number density of electrons in the pulsar wind follows directly 
from equation (|T3p and is given by 

"^^^ -r2iVe(r)5[7-7w(r)] 



d7df2dr 

rjL 



5[7-7w(r)] (26) 



47r7„(0)meC^ 
where we have assumed /^^ ~ 1- 

The radiative transport equation governing the specific intensity I{eout,r), 
including both emission by inverse Compton scattering and absorption by 
photon-photon pair production, is 

d-?^ 2 rdr 

— = eantrrieC T^ic + /— . (27) 
dr dr 

The quantity rjic is the inverse Compton emissivity, the number of photons 
emitted per second per unit volume per unit energy interval around eout'^eC^, 
which is given by 

1 r dNe dN^ 
Vic = ^ / d7 



mec"^ J d7dfidr deoutdt 

(28) 



1] dN^ 



An J 7„(0)(mec2)2cdeoutdt 

where dA^^/deoutdt is the differential rate of emission of photons by a single 
electron. 

The function Hideout, cos ip) which appears in the photon production rate, 
Eq. (|17D, is related to the quantity diV^/dtdeout[deindx'] derived by Ho & 
Epstein [1989, Eq. A14], and is defined by 

3^ 1 

?i(6out, COS^) = ^ [(i_;^_^)2 + ^2(l_^,2)]l/2 >< 

i + yJ + (^ + ^] + 

e (a2- 62)1/2 ^ ^ ey 



(a2- 62)1/2 J (a2_&2)3/2 



(29) 



with 

^ = 1 - P^x', A = eout/7w, £ = ein/7w, (30) 



24 



^, ^ /3w - cosV' 
1 — /9w cos -0 

and where 

a = e + £(l-2/o), 6 =-£5 (32) 



^ (i-a-0(i-a^pO .... 

^° (l-A-0' + /52(l-a;'2) ' ^ ^ 



^ /3(l-x-^)V^r + 2A(l-p)e-(l-pe)T/^ .^4. 

(1- A-e)2 + /32(l-a;'2) ' ^"^ 

and finally 

p = ^ . (35) 

f in 

However, this function gives the photon production rate only when t^^^ < 
tout < eiit where 

_ 7£(l-/jcOS^) _ 7£(l + £ + g,) 

(l + £ + g,) ' ^""^ " (e + 2£) ^^^^ 

with 

= {(3^ + + 2^8 00811^). (37) 

At these hmits the term + 2A(1 — p)^ — (1 — pO^]^''^ i'^ ^^e quantity 6 
becomes zero, and outside this range the photon production rate is zero. 

The differential optical depth due to pair production was given by Kirk et 
al. [1999, Eq. B4] In our notation it is 



if a; > 1 



, ^(1 -C') [(3-C^)lnfi4) -2C(2-C', ..^^^ 
otherwise 



where 



C^^^. (39) 
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